Dynamical mean-field theory for the anisotropic Kondo semiconductor: 
Temperature and magnetic field dependence 



(N 

o 

(N 

Oh- 
< 
(N 



-f— > 

CO 

c3 



i 

c 

o 
o 



> 
in 

O 



Takemi Yamada* and Yoshiaki Ono 
Department of Physics, Niigata University, Ikarashi, Nishi-ku, Niigata, 950-2181, Japan 

(Dated: March 2, 2013) 

We investigate the periodic Anderson model with fc-dependent of mixing reproducing the point 
nodes ol the hybridization gap by using the dynamical mean-field theory combined with the exact 
diagonalization method. At low temperature below a coherence temperature To, the imaginary 
part of the self-energy is found to be proportional to T 2 and the pseudogap with two characteristic 
energies Ai and A2 is clearly observed for T <C To, while the pseudogap is smeared with increasing 
T and then disappears at high temperature T > To due to the evolution of the imaginary self-energy. 
When the Coulomb interaction between / electrons U increases, A\, A2, and To together with T max 
at which the magnetic susceptibility is maximum decrease in proportion to the renormalization factor 
Z resulting in a heavy- fermion semiconductor with a large mass enhancement m* /m = Z~ x for large 
U. We also examine the effect of the external magnetic field H and find that the magnetization 
M shows two metamagnetic anomalies Hi and H2 corresponding to A\ and A2 which are reduced 
due to the effect of H together with Z. Remarkably, Z~ l is found to be largely enhanced due 
to H especially for Hi < H < H2, where the field induced heavy- fermion state is realized. The 
obtained results seem to be consistent with the experimental results observed in the anisotropic 
Kondo semiconductors such as CeNiSn. 



I. INTRODUCTION 

In the heavy- fermion systems, / electrons hybridize 
with conduction (c) electrons via the c-f mixing to form 
coherent quasiparticles with large effective mass, which 
is due to the effect of Coulomb interaction between / 
electrons at low temperature, while at high tempera- 
ture, / electrons are almost localized and scatter c elec- 
trons resulting in the Kondo effect. 1 The systems show 
various types of ground states including the so-called 
Kondo semiconductor, which exhibits an insulating be- 
havior at low temperature with highly reduced energy 
gap. With increasing temperature, the energy gap tends 
to be smeared, and then the system shows the behavior 
of incoherent metal at high temperature. Typical exam- 
ples of the Kondo semiconductors are SmBg, 2 YbBi2 3 
and Ce3Bi4Pt3, which have cubic crystal structures and 
possess well defined energy gaps of the orders of 100 K 
as observed in the measurements of the thermodynamic 
and transport properties. 

Another class of the Kondo semiconductors such as 
CeNiSn and its isostructural compounds with the or- 
thorhombic structure shows the behavior of anisotropic 
semiconductor or semimetal. 5-14 The longitudinal NMR 
relaxation rate 1 /T\ 10-12 and the Sommerfeld coefficient 
yr,i3,i4 are SU pp resS ed below 10 K indicating the devel- 
opment of a pseudogap in the density of states (DOS) 
at low temperature. Such pseudogap behavior is well 
accounted for by the V-shaped gap model with a resid- 
ual DOS 10 or the semimetallic model with nodes in the 
gap. 15 The V-shaped pseudogap in the DOS was directly 
observed below 10 K in the tunneling spectroscopy. 16 
As for the transport properties, the resistivity along the 
a-axis p a decreases with decreasing temperature as ex- 
pected by the semimetallic model, 15 while pb and p c 
slightly increase below 3 K. 7 The inelastic neutron scat- 



tering experiments revealed the existence of anisotropic 
magnetic excitations. 17-19 Anisotropic pseudogap prop- 
erties were also observed in the magnetization and the 
magnetoresistance. 20-22 Despite the intense efforts, it is 
still controversial whether CeNiSn has a zero DOS just at 
the Fermi level or a dip structure DOS around the Fermi 
level within the experimental uncertainties. 

A remarkable feature of the anisotropic Kondo semi- 
conductors is the significant temperature dependence of 
the pseudogap. In the tunneling 16 and the photoemis- 
sion spectroscopies, 23 ' 24 the pseudogap in the DOS is 
clearly observed at low temperature, while it vanishes at 
high temperature. Such temperature dependent pseudo- 
gap cannot be explained with a simple rigid-band model. 
Therefore the electron correlation effect is considered to 
be crucial for the temperature dependence of the pseudo- 
gap together with the large reduction of the gap width. 
More recently, possible long-range ordered states in the 
anisotropic Kondo semiconductors have also been exten- 
sively investigated with the effects of the pressure 25 ' 26 
and the doping. 27-29 

Many theoretical studies for the Kondo semiconductor 
have been made on the basis of the periodic Anderson 
model (PAM) 30-32 with fc-independent c-f mixing repro- 
ducing the isotropic hybridization gap by means of vari- 
ous methods such as the Gutzwiller approximation, 33 ' 34 
the slave-boson mean-field theory, 35-37 the noncrossing 
approximation, 38 the 1/N expansion 39-41 and the dy- 
namical mean-field theory (DMFT). 42-47 These studies 
have shown that, due to the strong correlation effect, the 
hybridization gap is highly reduced to form a renormal- 
ized gap 33-37 which is clearly observed at low tempera- 
ture but disappears at high temperature. 38-45 A magnetic 
field induced insulator to metal transition has also been 
observed at low temperature. 46,47 The obtained results 
are consistent with the experimental results observed in 
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the isotropic Kondo semiconductors such as Ce3Bi 4 Pt 3 
and YbB 12 . 

As for the anisotropic Kondo semiconductors such as 
CcNiSn, the fe-dependence of the c-f mixing is consid- 
ered to be important in addition to the strong corre- 
lation effect. The fc-dependent c-f mixing originates 
from the crystal electric field (CEF) ground states of 
/ electrons 31 ' 32,48 and yields the specific DOS with 
the pseudogap structure. 15,49-51 The PAM with the fc- 
dependent c-f mixing has been studied by using the 
Gutzwiller approximation 15 and the slave-boson mean- 
field theory, which reproduce the highly reduced pseu- 
dogap and well explain the thermodynamic and transport 
properties of the anisotropic Kondo semiconductors at 
low temperature. However, the temperature dependence 
of the pseudogap, which is directly observed in the tun- 
neling and photoemission spectroscopies, 16,23,24 together 
with the magnetic field dependence was not discussed 
there. 

The purpose of this paper is to elucidate the effects of 
temperature and magnetic field on the electronic states 
of the anisotropic Kondo semiconductors. For this pur- 
pose, we study the PAM with fc-dependent c-f mixing 
by using the DMFT which becomes exact in the limit 
of infinite spatial dimensions and is expected to be a 
good approximation in three dimensions. The DMFT 
is known to describe well the strongly correlated electron 
systems over the whole parameter regime of temperature, 
magnetic field and frequency, and has been extensively 
developed for the PAM with fc-independent c-f mixing 
to describe the heavy-fermion systems and the Kondo 
semiconductors. 42-47 In the previous work, we have em- 
ployed the DMFT combined with the exact diagonaliza- 
tion (ED) method for the PAM with the fc-dependent c-f 
mixing and have obtained the magnetic field dependence 
of the electronic state which well accounts for the meta- 
magnetic behavior observed in CeRu2Si2- 52 The present 
paper is a straight forward extension of the previous work 
for the case with the anisotropic Kondo semiconductors 
such as CcNiSn. 

In this paper, we investigate the anisotropic Kondo 
semiconductor on the basis of the PAM with the 
fc-dependent c-f mixing at half-filling by using the 
DMFT+ED method. 52 The physical quantities are cal- 
culated systematically over the wide parameter regime of 
temperature T, magnetic field H and Coulomb interac- 
tion U between / electrons. The paper is organized as 
follows: in Sec. II, we present the Hamiltonian of the 
PAM with the fc-dependent c-f mixing and the formula- 
tion of the DMFT+ED method. In Sec. Ill, we show the 
results of the physical quantities for H = 0, the renor- 
malized DOS, the magnetic and charge susceptibilities, 
the renormalization factor and the imaginary part of the 
self-energy as functions of U and T. In Sec. IV, we 
present the H dependence of the physical quantities, the 
magnetization and the renormalization factor for various 
U and T. In Sec. V, the paper is ended with a summary 
together with discussions where the present results are 



compared with the previous theoretical results and the 
experimental results. 

II. MODEL AND FORMULATION 

A. Model Hamiltonian 

Our model Hamiltonian of the PAM with fc-dependent 
c-f mixing 15,48,51 consists of the conduction electron 
term H c , the / electron term Hf and the c-f mixing 
term H c f as follows 

H = H c + Hf + H c f , (1) 

H c = ^2e k cl a c k a, (2) 

ha 

H f = Y e /™ n im + U Y n i+ n L > ( 3 ) 

im i 
H cf = ( V kmaft m Ck* + H.C.) , (4) 

where c km is a creation operator for a c electron with the 
wave vector fc and the spin a =t, 4-, f\ m is that for a / 
electron with the lowest Kramers doublet state m = ± at 
site i and n{ m = f\ m fi m - e k (e fm ) is the energy for the c 
(/) electron, V kma is the c-f mixing matrix element and 
U is the Coulomb interaction between / electrons. In Eq. 
(4), the effect of the external magnetic field H is included 
only in the / electrons as e/ m = e/ — mH, because the g 
value for the / electron is known to be much larger than 
that for the c electron. 

In this paper, we assume that the lowest Kramers dou- 
blet state under the CEF is J z = ±3/2, which is referred 
to as m = ± and the c electron state is simply given by 
the plane wave. Such a model was originally developed 
by Ikeda and Miyake to describe the electronic state of 
the anisotropic Kondo semiconductors such as CeNiSn. 15 

In this case, the c-f mixing matrix element is given by 
[see Eq. (A10) in Appendix] 

f v k+t V k+i \_ ( -aY 31 (n k ) bY 32 (n k ) \ 

\ v fc - t v k ^ ) - Vc f \ -bY 3 _ 2 (n k ) ay 3 -i(n fc ) / ' 

with a = \J^- and b = \J~^-, where Y pq (Q k ) is a spher- 
ical harmonics with the argument of the solid angle Q k 
of the wave vector fc, and 14/ is the c-f mixing strength 
defined in Eq. (A6) in Appendix, which is a parameter 
in our model. 

B. c-f hybridized bands for U = 

In the nonintcracting case with U = 0, the Hamilto- 
nian Eqs. (1) — (4) with Eq. (5) is diagonalized to yield 
the c-f hybridized bands with the energies, 

E ( k 2 = \ Urn + e k ± V /(e /m -e fc ) 2 +4/ fe ") (6) 
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FIG. 1. The noninteracting DOS for / and c electrons p/(e) 
and p c (e) near e = for U = 0, e/ m = 0, V c / = 0.5 and 
D = 2. 



with /fe = |Vfc m(T | 2 , where depends only on the 
z component of the unit fc-vector, k z — k z /\k\, and is 
explicitly given by (see Eq. (A13) in the Appendix) 



/„ = -^(1-^(1 + 15^ 



(7) 



Then, the hybridization gap between the upper and the 
lower hybridized bands has nodes on the k z axis with 
k z = ±1 resulting in a pseudogap structure of the DOS 
as shown in Fig. 1, where two characteristic energies of 
the pseudogap A\ and A 2 originates from the minimum 

of Ik at k z = and the maximum of If, at k z = ± ^J~^, 
respectively. This pseudogap is found to well reproduce 
the anisotropic Kondo semiconductor CcNiSn, 15 where 
the characteristic temperature dependence of the specific 
heat and the NMR relaxation rate at low temperature is 
well accounted for by the pseudogap. We note that the 
/ electrons DOS at the Fermi level is finite as shown in 
Fig. 1, and then, the resistivity shows a metallic behav- 
ior at low temperature as observed in CeNiSn. 15 This is a 
striking contrast to the case with the fc-independent c-f 
mixing that yields a finite hybridization gap reproducing 
the isotropic Kondo semiconductor, where the resistiv- 
ity shows a semiconducting behavior at low temperature. 
Here and hereafter, we assume the bare c-DOS to be a 
rectangular DOS with the band width 2D centered at 
e = 0. 



C. DMFT+ED formalism 

In the DMFT, 45-47,52 the lattice model is mapped onto 
an effective impurity model embedded in an effective 
medium which is to be determined self-consistently. To 
solve the effective impurity model, we employ the ED 
method for a finite-size cluster given by the following 



Hamiltonian, 



-ffimp — 



N.-l 



-1 m 



N.-l 



2 (Vi m a\ m ai- lm + H.c.) , (8) 



1 = 1 rn 



where a\ m is a creation operator for an electron with m = 
± for the impurity site I = and that for the effective 
medium sites I = 1, • • • , iV s — 1, respectively, and ni m — 
a\ m aim- U is the Coulomb interaction between electrons 
on the impurity site and is set to be the same value of 
U in the original lattice Hamiltonian Eq. (3). A set of 
parameters {e; m , Vim} is so-called Weiss field parameters 
(WFPs) , which represents the effective medium and is to 
be determined self-consistently. 

In the noninteracting case with U = 0, the impurity 
Green's function is written with the WFPs as 



\Vlr, 



1=1 



z v - eir. 



(9) 



where z v — i{2v + 1)ttT is the Matsubara frequency 
with the temperature T. For finite U, we solve the iV s - 
site Hamiltonian (8) by using the Householder ED algo- 
rithm to obtain the impurity Green's function G m (z l/ ) = 
(Gmizu)^ 1 — S m (z !/ )) _1 together with the impurity self- 
energy X m (z„). Then, the self-consistency condition, 
where the impurity Green's function coincides with the 
local / electron Green's function G^z^) of the original 
PAM in Eq. (1) with the same self-energy Y, m {z v ), is 
given by 



1 



gO l (^)- 1 -S m (z,)' 



(10) 



Substituting Eq. (7) into G^Zy) in Eq. (10), we obtain 
a more explicit expression for G' m {z) as 



G f m (z) 



Z e fm ~ S m (z) 



x — / dk z 7 fc ln 



(z - ef m - E m (z)) 2 
C + D 



C-D 



(11) 



with C = z v - I k /(z u - ef m - T, m (z)). 

In the explicit calculation to obtain the DMFT+ED 
solution, the following procedures are carried out: (i) 
under given WFPs {e°^, V^ d }, the AT s -site Hamiltonian 
Eq. (8) is solved by using the Householder ED algorithm 
to obtain eigen values and eigen vectors, (ii) From the 
eigen values and eigen vectors, the self-energy T lm {z u ) 
is calculated, (hi) Substituting T, m (z u ) into Eq. (10), 
new WFPs ief m w , V,^ w } are determined so as to satisfy 
the self-consistency condition Eq. (10) with Eq. (9) as 
possible. The steps (i) — (hi) are iterated until the old and 
new WFPs coincide with each other. 
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III. RESULTS FOR H = 

In this section, we show the U and T dependence of 
physical quantities in the absence of the external mag- 
netic field, H = 0. We restrict ourselves only to the 
half-filling case with the particle-hole symmetry, where 
we set 6f = — U/2 and the c and / electron numbers per 
site are given by (n c ) = (n?) = 1. In addition, we as- 
sume that the system is in the paramagnetic state with 
{n + ) = (n_) = ^, and then the physical quantities are 
independent of m. Here and hereafter, the parameters 
are set to as follows: the half c-band width D = 2, the 
c-f mixing strength V c j = 0.5, and the site number of 
the effective impurity model N s — 6. 53 

A. Renormalized /-DOS 

Figure 2(a) shows the renormalized /-DOS pf(e) near 
the Fermi level e = for several values of U at a low 
temperature T = 0.001. In this calculation, we perform 
the analytic continuation of the self-energy S m (z) from 
the imaginary frequency to the real frequency by using 
the Pade approximation, and then substitute E m (e+i0+) 
into Eq. (11) to obtain //(e) = p/„(e) = — ilm G^e + 
i0_|_). When U increases, the renormalized pseudogap 
with the renormalized characteristic energies Ax and A 2 , 
which correspond to A\ and A2 for U = shown in Fig. 
1, decreases together with decrease in the quasiparticle 
band width resulting in the highly reduced pseudogap ac- 
companied by the heavy-fermion bands for large U. We 
also find that the spectral weight due to the quasiparticle 
bands decreases with decreasing the quasiparticle band 
width, while that due to the broad peaks correspond- 
ing to the Hubbard- like bands around e/ = — U/2 and 
ef + U = U/2 increases (not shown). We note that, as 
the self-energy E m (z) is independent of k in the DMFT, 
pf (0) at T = is unchanged by U and has a finite value 
(see also Fig. 1) resulting in a metallic behavior at low 
temperature as mentioned in Sec. II B. 

One of the most remarkable features of the Kondo 
semiconductor is that the gap structure largely depends 
on temperature in contrast to the case with the ordi- 
nary semiconductor. In Fig. 2(b), we plot the renor- 
malized /-DOS p f (e) near the Fermi level e = for sev- 
eral values of T at U = 1.8. At low temperature, we 
observe a clear pseudogap structure with the renormal- 
ized characteristic energies A\ and A 2 . With increas- 
ing T, the pseudogap structure is found to be smeared, 
and then finally disappears at high temperature above 
the so-called coherence temperature To, where the T de- 
pendence of pf(e) is mainly caused by the evolution of 
the imaginary part of the self-energy ImS(e) that be- 
comes large for T > To as explicitly shown in Sec. IIIC. 
Such T dependence of p*(e) has been observed in the 
tunneling 16 and photoemission 23 ' 24 spectroscopies for the 
Kondo semiconductors. 

To see the renormalized characteristic energies A\ and 



Ai more explicitly, we plot the energy derivative of the 
/-DOS dp f (e)/de for several values of U at T = 0.001 as 
shown in Fig. 3. The two significant peaks of dpf(e)/de 
corresponding to A±/2 and A2/2 are clearly observed and 
found to decrease with increasing U. In Fig. 4, the U 
dependence of A\ and A 2 are plotted together with the 
renormalization factor Z calculated from the real part of 

the self-energy as Z = Z m = II — ^Re£ m (e) ) 

V € — 0/ 

When U increases, A\ and A2 decrease with decreasing 
Z resulting in the highly reduced pseudogap. 

As mentioned in Sec. II B, the pseudogap energies A\ 
and A2 are given by the band energies with k z = and 
k z = ±^/7/15, respectively. Then, A\ and A 2 are ex- 
pected to be given by the corresponding energies of the 
renormalized quasiparticle bands whose widths are re- 
duced by the renormalization factor Z. In fact, A\ and 
A2 arc found to be in good agreement with A\Z and 
A 2 Z, respectively, and are highly reduced in proportion 
to Z <C 1 for large U as shown in Fig. 4. The highly 
reduced pseudogap observed for large U is accompanied 
by the heavy fermions with the large mass enhancement 
factor m*/m = Z^ 1 ^> 1 resulting in the heavy-fermion 
semiconductor with the nodal gap structure. 

B. Magnetic and charge susceptibilities 

In Fig. 5, we show the uniform and the local compo- 
nents of the magnetic susceptibilities for / electrons x™ 1 
and x m c as functions of T for U = 0, 1 and 1.5, where x™ 1 
is calculated from the magnetization M = (nl) — (n{_) in 
the presence of a small external magnetic field H = 0.01 
as xUf 1 = M/H and Xm° 1S calculated from the eigenval- 
ues and the eigenvectors in the effective impurity model 
for H = using the standard linear-response formula- 
tion. We can see both of the magnetic susceptibilities 
agree well with each other in the case with the present 
model where the specific nesting vector responsible for 
the strong fe-dependencc of the magnetic susceptibility is 
absent. 

For U = 0, x™ 1 (together with Xm C ) shows a maximum 
at a certain temperature T max which roughly corresponds 
to Ai/2 and is described by the Pauli paramagnetism 
with the pseudogap structure of the /-DOS (see Fig. 1). 
As shown in Fig. 5, x™ 1 (Xm C ) i s enhanced for U = 1 
and U — 1.5 due to the correlation effect resulting in 
an enhanced Pauli paramagnetism at low temperature 
below the coherence temperature T , where the pseudo- 
gap structure of the renormalized /-DOS is observed [see 
Fig. 2(b)], and then, x™ 1 (Xm°) shows a maximum at 
Tmax ~ ^i/2 due to the pseudogap structure. On the 
other hand, at high temperature T Z To, xUf' (x m c ) ex_ 
hibits a Curie-law behavior x™ 1 ~ Xm° ~ 1 /T where the 
/ electron is considered to be almost localized and then 
the pseudogap structure in the renormalized /-DOS is 
found to disappear as shown in Fig. 2(b). 

To see the correlation effect due to U systematically, 



5 



fa) 



(a) [/dependence at 7=0.001 



-0.1 



fa 




U=2.0 
U=l.9 
£/=1.8 
£/=1.7 
£7=1.6 
£7=1.5 
£7=1.4 
£7=1.3 
£7=1.2 
£7=1.1 
£7=1.0 



-0.05 



7=0.040 

7=0.017 
r=o.oi5 

7=0.013 

r=o.on 
7=0.009 

7=0.007 

7=0.005 
7=0.003 
7=0.001 



FIG. 2. (Color online) The renormalized /-DOS fr(e) near 
the Fermi level e — for several values of [7 at T = 0.001 (a) 
and for several values of T at £/ = 1.8. 



we plot the T dependence of Xm C f° r several values of £7 in 
Fig. 6. When £7 increases, x|^ c increases to show the en- 
hanced Pauli paramagnetism at low temperature T ;$ To 
together with the Curie-law behavior at high temperature 
T Z To, where A\ decreases with increasing £7 as shown 
in Fig. 4. In Fig. 6, we also plot the T dependence of the 
local charge susceptibility for / electrons Xc°°j which is 
also calculated from the eigenvalues and the eigenvectors 
in the effective impurity model using the standard linear- 
response formulation. x|° c monotonically decreases with 
increasing £7 and is largely suppressed for large £7 corre- 
sponding to the Kondo regime. 



C. Quasiparticle lifetime 
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FIG. 3. (Color online) The energy derivative of /-DOS 
dp '(e) /de for several values of £7. 



that — Im£(i0+) obeys the T 2 dependence at low tem- 
perature below the coherence temperature T$ and is well 
reproduced by a fitting function cq + c (T/Tq) , where cq 
the value at T — expected to be zero in the Fermi liquid 
theory and c is the value at T = lo- 
in the present numerical calculation, Co is small (co <C 
c) but finite due to the effect of finite size N s , and is found 
to decrease with increasing N s as approaching Co — > 
for N s — > oo. Then, the inverse lifetime is estimated 
as l/2r = — Im£(i0 + ) — c and is plotted as a function 
of T 2 for several values of £7 in Fig. 7. We can see 
that l/2r thus obtained is in good agreement with the 
fitting function l/2r = c(T/T Q ) 2 (dotted lines), where 
we set c = 0.04 above which l/2r is found to deviate 
from the T 2 dependence and we determine To so as to fit 
l/2r to the fitting function as possible for each £7. When 
£7 increases, the coherence temperature To decreases in 
proportion to Z as explicitly shown in the next section. 



The quasiparticle lifetime r is known to be related 
to the imaginary part of the self-energy as l/2r = 
— Iml!(i0-|_). To obtain ImS(i0+), we perform the ana- 
lytic continuation of the self-energy S m (z) from the imag- 
inary frequency to the real frequency by using the Pade 
approximation as mentioned in Sec. Ill A. It is found 



D. Characteristic temperatures 

In Fig. 8, we plot the characteristic temperatures, 
Tnax obtained in Sec. Ill B and T obtained in Sec. Ill C, 
as functions of £7 together with the characteristic ener- 
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FIG. 4. (Color online) [/-dependence of the characteristic 
energies of the renormalized pseudogap A\ and A2 and the 
renormalization factor Z together with A\Z and A\Z, where 
A\ and A2 are the bare characteristic energies of the pseudo- 
gap. 
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FIG. 6. (Color online) T dependence of the local magnetic 



and charge susceptibilities Xm C (T) and 
for several values of U. 



Xc(T), respectively, 




FIG. 5. (Color online) T dependence of the uniform and the 
local components of the magnetic susceptibilities for / elec- 
trons Xm 1 and Xm c for U = and U = 1. 



gies of the renormalized pseudogap A\ and Z\ 2 obtained 
in Sec. Ill A. We can see that To monotonically decreases 
with increasing U and is roughly proportional to Z for 
large U, where To ~ A± ~ A\Z (see also Fig. 4). Then, 
the T 2 coefficient of the inverse lifetime A — c/Tq is pro- 
portional to Z~ 2 — (m*/m) 2 for large U as expected 
from the Fermi liquid theory. We note that, for small 
U, the inverse lifetime, i. e., the imaginary part of the 
self-energy, can be obtained from the second-order per- 
turbation with respect to U and is proportional to U 2 , 
and then T oc U' 1 as observed in Fig. 8. 

At low temperature T < To, the imaginary part of the 
self-energy is sufficiently small to obtain the well defined 
quasiparticles, which yield the pseudogap structure of the 
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FIG. 7. (Color online) T 2 -dependence of the inverse lifetime 
for several values of U. The dotted lines represent the fitting 
functions l/2r = c(T/T ) 2 (dotted lines) with c = 0.04 and 
the coherence temperature To which depends on U. 



renormalized /-DOS as shown in Fig. 2 (b). Therefore 
the T dependence of the physical quantities such as Xm 
is well described by the quasiparticle band for T < To, 
where Xm shows a maximum at T max ~ A\/2, which 
is smaller than To, due to the pseudogap structure of 
the renormalized /-DOS as shown in Fig. 5. When U 
increases, T max ~ A\/2 decreases in proportion to Z as 
A x - AxZ (see Fig. 4). 
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FIG. 8. (Color online) U dependence of the characteristic 
temperatures T max and To together with the characteristic 
energies of the renormalized pseudogap Ai and A-z- 



IV. RESULTS FOR H / 

In this section, we examine the effect of the external 
magnetic field H especially focused on the metamagnetic 
behavior. Due to the particle-hole symmetry with £/ = 
—{7/2, the self-energy E m (z) is independent of m even 
for H ^ 0, and then, Z + = Z_ = Z. 



A. Magnetization and differential susceptibility 

Figure. 9 shows the H dependence of the magnetiza- 
tion M = (n+) — (n/) and that of the differential sus- 
ceptibility dM/dH for several values of U and T. At low 
temperature, two metamagnetic anomalies in M are ob- 
served at critical magnetic fields H = Hi and H 2 as 
shown in Figs. 9(a) — (c), and the corresponding two 
sharp peaks in dM/dH are observed at H = Hi and 
H2 as shown in Figs. 9(d) — (f). We find that the critical 
magnetic fields Hi and Hi are given by the following rela- 
tions: Hi = Ax{Hx)/2 and H 2 = A 2 (H 2 )/2, respectively, 
where Ai(H) and A 2 (H) are the renormalized pseudo- 
gap energies whose values for H = are shown in Fig. 4 
and largely depend on the external magnetic field H as 
explicitly shown in the next section. 

When U increases, Ai(H) and A 2 (H) decrease for 
a given value of H as explicitly shown in Fig. 4 for 
H = 0. Therefore Hi and H 2 decrease with increasing 
U as shown in Figs. 9(d) — (f). When T increases, the 
two metamagnetic peaks in dM/dH clearly observed at 
T < Tmax are smeared, and then, merge into a broad peak 
at T max ~ T < T , where T max and T are shown in Fig. 
8. The metamagnetic behavior finally disappears to show 
the monotonic M-H curve similar to the case with the 
localized spins at high temperature T > T , where the 
pseudogap structure disappears as shown in Fig. 2(b) 




FIG. 9. (Color online) H dependence of the magnetization 
M at U = 1.6 (a), U = 1.8 (b), and U = 2.0 (c) and that of 
the differential susceptibility dM/dH a,tU = 1.6 (d), U — 1.8 
(e), and U = 2.0 (f) for several values of T. 



and the magnetic susceptibility shows the Curie-law be- 
havior as shown in Fig. 6. 



B. Inverse renormalization factor 

As mentioned in the previous subsection, the renor- 
malized pseudogap energies Ai(H) and A 2 {H) largely 
depend on the external magnetic field H at low temper- 
ature. This is mainly caused by the H dependence of the 
renormalization factor Z(H) as Ai{H) « Z(H)Ai and 
A 2 (H) fa Z(H)A 2 . In Fig. 10, we plot the H dependence 
of the inverse renormalization factor Z~ x for several val- 
ues of U and T. At low temperature, Z~ x increases with 
increasing H for H ^ Hi, where it shows an abrupt in- 
crease at H — Hi, and shows a peak for Hi < H < H 2 , 
and then decreases for H > H 2 , where it shows a kink 
at H = H 2 . This H dependence is caused by the pseu- 
dogap structure of the renormalized /-DOS p?(e) at low 
temperature (see Fig. 2 for H = 0), where the /-DOS 
at the chemical potential p^(0) is small due to the pseu- 
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dogap for H < Hi, while it shows an abrupt increase 
at H ~ Hi and becomes large for H\ < H < iJ 2 , and 
then gradually decreases with increasing H for iJ > ii^- 
Therefore the correlation effect between the / electrons 
is enhanced due to H especially for Hi < H < H 2 re- 
sulting in the enhancement of m*/m = Z . Then, we 
observe a remarkable field induced heavy-fermion state 
forHi<H<H 2 . 

When T increases, Z" 1 decreases for Hi < H < H 2 , 
while it increases for H ^ Hi and H Z H 2 . Therefore, the 
peak structure of Z~ x with pronounced anomalies at Hi 
and H 2 clearly observed at low temperature T < T max 
is smeared to show a broad peak around Hi < H < 
H 2 at T max < T < T , and then finally disappears to 
show the monotonically decreasing function of H at high 
temperature T > To, as similar to the case with dM/dH 
[see Figs. 9(d) -(f)]. 

Finally, we plot the differential susceptibility dM/dH 
and the inverse renormalization factor Z~ 1 as functions 
of H for various U at a low temperature T = 0.001 in 
Figs. 11(a) and 11(b). At H = Hi, dM/dH shows a 
sharp peak and Z~ x shows an abrupt increase, while, at 
H = H 2 , dM/dH shows a cusp and Z' 1 shows a kink. 
For Hi < H < H 2 , both of dM/dH and Z" 1 are largely 
enhanced as compared to those values for H = 0. Thus, 
we find that the pseudogap structure due to the point 
nodes of the hybridization gap in the PAM is responsi- 
ble for the remarkable field induced heavy-fermion state 
accompanied by the metamagnetism. 



V. SUMMARY AND DISCUSSIONS 

In summary, we have investigated the PAM with the 
fc-dependent c-f mixing that reproduces the point nodes 
of the hybridization gap simulating the pseudogap struc- 
ture of the anisotropic Kondo semiconductors by using 
the DMFT+ED method. The physical quantities have 
been calculated systematically over the entire range of 
the parameters: T, U, and H. What we have found are 
as follows: (1) at low temperature below the coherence 
temperature To, the imaginary part of the self-energy is 
proportional to T 2 , where the pseudogap with two char- 
acteristic energies Ai and A 2 is observed. The magnetic 
susceptibility shows the enhanced Pauli paramagnetic be- 
havior with a maximum at T max (< To) due to the effect 
of the pseudogap. When U increases, Ai, A 2 , T and 
T max decrease in proportion to the renormalization fac- 
tor Z resulting in a heavy-fermion semiconductor with 
a large mass enhancement m*/m = Z^ 1 for large U. 
(2) In the presence of the external magnetic field H at 
low temperature T < To, the magnetization M shows 
two metamagnetic anomalies Hi and H 2 corresponding 
to A i and A 2 that are reduced due to the effect of H 
together with Z. Remarkably, Z~ x is largely enhanced 
due to H especially for Hi < H < H 2 , where the field 
induced heavy-fermion state is realized. (3) When T in- 
creases, the pseudogap together with the metamagnetic 
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FIG. 10. (Color online) H dependence of the inverse renor- 
malization factor Z" 1 for several values of T at U = 1.6 (a), 
U = 1.8 (b) and U = 2.0 (c). 



anomalies is smeared due to the evolution of the imagi- 
nary self-energy and finally disappears at high tempera- 
ture T <; To, where the magnetic susceptibility shows the 
Curie-law behavior. 

The present DMFT results are consistent with the pre- 
vious results from the Gutzwiller approximation 15 and 
the slave-boson mean-field theory 49 at low temperature 
T <C To where the renormalized pseudogap is clearly ob- 
served. On the other hand, with increasing T, the de- 
viation between the present and the previous results in- 
creases appreciably due to the evolution of the imaginary 
self-energy resulting in the smearing of the pseudogap, 
and then, the discrepancy becomes significant at high 
temperature T ^ To, where the pseudogap disappears in 
the present study in contrast to the previous studies with 
the rigid pseudogap structure. 

In the previous DMFT studies, the PAM has been ex- 
tensively investigated in the case with the fc-independent 
c-f mixing, 45-47 where the fully opened hybridization gap 
A is found to be highly reduced to a renormalized value 
A w ZA due to the strong correlation effect. At the half- 
filling with H = 0, the Fermi level sits in the renormal- 
ized hybridization gap resulting in the insulating ground 
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FIG. 11. (Color online) H dependence of the differential sus- 
ceptibility dM/dH (a) and inverse renormalization factor Z 
(b) for several values of U at T = 0.001. 



state, i.e., the isotropic Kondo semiconductor. When 
H increases, the DOS is split due to the Zeeman split- 
ting, and then, the Fermi level enters the upper (lower) 
hybridized band of the majority (minority) spin above 
a critical magnetic field H c ~ A/2 at which the field 
induced insulator-metal transition takes place. 46 Corre- 
spondingly, Z~ 4 is found to be almost independent of H 
below H c , while Z~ x shows a discontinuous increase at 
H = H c , 46 where the DOS at the Fermi level changes 
from zero to a finite value and then the correlation effect 
is enhanced. For H > H c , Z^ 1 monotonically decreases 
with increasing H as the Kondo effect due to the local 
spin fluctuation is suppressed by the magnetic field. At 
finite temperature, the discontinuous increase in Z^ 1 ob- 
served at T — changes into a broad peak around H c 
due to the thermal broadening effect. 47 

In the present PAM with the fc-dependent c-f mixing 
reproducing the anisotropic Kondo semiconductor, the 
Fermi level sits in the dip of the DOS at the half-filling in- 
stead of in the fully opened hybridization gap with the fc- 
independent c-f mixing. Therefore the DOS at the Fermi 
level is finite even for H — and increases with increasing 



H toward H\. Correspondingly, the mass enhancement 
factor Z~ x gradually increases with increasing H toward 
Hi, where the correlation effect is enhanced due to the 
large DOS at the Fermi level. This is the origin of the 
field induced heavy-fermion state in the present model. 
For H ^> Hi, Z^ 1 gradually decreases with increasing H 
as the Kondo effect is suppressed by the magnetic field. 
We note that, in the present study, we concentrate only 
on the particle-hole symmetric case, where Z m is inde- 
pendent of to: Z + = Z- = Z, even for finite H . In the 
case without the particle-hole symmetry, the m depen- 
dence of Z m for finite H is considered to be significant. 
In fact, such to dependence of Z m was obtained in our 
previous DMFT study for the PAM away from the half- 
filling reproducing the heavy-fermion metamagnetism. 52 
Finally, we compare our theoretical results with the 
experimental results in the anisotropic Kondo semicon- 
ductors such as CeNiSn. At low temperature T < Tq, 
the pseudogap structure due to the fc-dependent c-f mix- 
ing well describe the characteristic T dependence of the 
physical quantities observed in CeNiSn as previously ob- 
tained by Ikeda and Miyake. 15 In addition, the present 
results account for the T dependence also at high tem- 
perature T ^ Tq, e.g., the Curie law behavior of the 
magnetic susceptibility. Most significantly, the T depen- 
dence of the /-DOS in the present results is consistent 
with the experimental results from the tunneling 16 and 
photoemission 23,24 spectroscopies, where the pseudogap 
structure is clearly observed at low temperature, while 
it disappears at high temperature. The metamagnetic 
anomalies at Hi and Hi are also consistent with the 
experimental observations. 20-22 More recently, the high 
magnetic field measurements exceeding 50 T (see Ref. 
22) have revealed that the existence of a third metamag- 
netic anomaly at H^(> H2) which has not been obtained 
in the present study. To discuss the anomaly at H3, the 
effect of the excited CEF levels, which are not included in 
the present study, is considered to be important. There- 
fore, we need further investigations to describe the T and 
H dependencies of the Kondo semiconductors by using 
more realistic models including the excited CEF levels 
together with the realistic band structure. 
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Appendix A: fc-dependent c-f mixing matrix 

In this appendix, we note the general expressions of the 
fc-dependent c-f mixing term. 15 ' 48 ' 50,51 Here c electron 
state is a direct product of the plane wave state and the 
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spin state x<r(£) as 



Vw(r, = <r£|fc<r) = —e^ Xa {t). (Al) 
V i v 



The plane wave expansion around a position vector r is 
given by 

e ik-r = ^J2i P j P (kr)Y p * q (e k , <j> k )Y pq {6, <f>), (A2) 

pq 

where j p {kr) is a spherical Bessel function and Y pq (9,4>) 
is a spherical harmonics with the argument of the solid 
angle £l r of the vector r or fl k of the wave vector k. 
Therefore, the c electron wave function is described by 

x ^^(fcr^F^^^y^^Jx^O, (A3) 



where r j is a vector from a site i?i , = r — Ri. 

The / electron states are characterized by the localized 
atomic orbital around a site Ri . The LS coupling in the 
Ce compounds splits the 14-folded wave functions to the 
excited J = 7/2 octet and the ground J = 5/2 sextet. 
In Ce compounds, the sextet in J = 5/2 are usually 
split into three Kramers doublets by CEF. Each Kramers 
doublet is expanded with a linear combination of J 
multiplets as below 

,7 

M r >0 = H\») = E (r£\JM)(JM\ri, 

M=-J 

= E b t*MaMqaRnp(r)Y pq (9,4>)Xa(0> 
Mqa 

where R np (r) is the radial function with quantum number 
(n,p) and auqa = (npqa\JM) is the Clcbsch-Gordan 
coefficient and 6 m m = (JM\/j,) is a coefficient with the 
linear combination of some of \JM). Consequently, / 
electron wave function around a site Ri is described by 

Vv( r »>0 = E b nMaMqaRnp(ri)Y pq (n ri )Xa(0 , 
Mqa 



where au q a = — &\J -^-j^-5 q M-% for the J = 5/2 sex- 
tet. By using the above wave functions, c-f mixing be- 
tween the states \ka) and can be directly calculated 
as below, 



(A4) 



where V kll(7 by 

V kflc7 = V&rVknp E h iiMa M qaY pq {n k ), (A5) 

Vfe„ p = ViTT(-i) p J dr, r 2 j p (kri)v(r t )R np (n), (A6) 

where the mixing strength Vknp = K,/ is a treated as a 
parameter of our model. Thus, the c-f mixing Hamil- 
tonian between the states \ka) and £f c / is written 

by 



H cf = E (Vfe^C^/fep + H.c. 



kfj,a 



(A7) 



In general, the J = 5/2 sextets are split into three 
Kramers doublets except for the Tg state in the cubic 
CEF. When we consider the only lowest Kramers dou- 
blets denoted by the pseudo spin states \x = ±, the c-f 
mixing matrix V k with their elements V kflrT scaled by V c f 
is a 2 x 2 matrix, which is written by 



V k 



J_(v k+t v k+l \ 



(A8) 



For example, when \x = ±5/2, 
V k = 



'^Y 32 (Q k ) 
'^Y 3 _ 3 (n k ) 



/2 4?Y 33 (n k ) 
^Y 3 _ 2 {n k ) 



(A9) 




(A10) 



(AH) 



The important quantity for the fe-dependent c-f mixing 
is I k as mentioned in the Sec. II B, and these are written 

by 

^ = ±5/2 J* = yV£(l-#) 2 , (A12) 

M = ±3/2 I k = ^V c 2 f (l-k 2 z )(l + 15k 2 z ), (A13) 

/i = ±1/2 I k = - A V c 2 f (l-2k 2 z + bkt). (A14) 
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